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Abstract 

Deterministic diffusion in temporally oscillating convection is studied 
for particles with finite mass. The particles are assumed to obey a simple 
dissipative dynamical system and the particle diffusion is induced by the 
strange attractor. The diffusion constants are numerically calculated for 
convection models with free and rigid boundary conditions. 
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The motion of a fluid particle in the velocity field v(x) is determined by the 
differential equation. 

f=v M (1) 

with the initial condition x(0). Even for regular velocity fields, the particle 
motion can become chaotic and it is called the Lagrangian chaos jj], §]. For an 
incompressible fluid, the dynamical system (1) is a conservative system. In two 
dimensions, the evolution equations become 

dx dip dz dip , x 
dt dz ' dt dx' 

where ip is the stream function and x = (x, z). Equation (2) is a Hamiltonian 
system. Solomon and Gollub studied experimentally and numerically the par- 
ticle diffusion in Rayleigh-Benard convection which oscillates temporally |3| . In 
the oscillating convection, the velocity field is derived from a stream function 
ip = ^4/fcsin{fca; + B sm(uit)}W (z) , where k is the wavenumber, A is the maxi- 
mum vertical velocity, and B and w are the amplitude and the frequency of the 
oscillation and W(z) is a function which satisfies the boundary conditions at 
the top and bottom surfaces (z = 1 and 0). 

In this report, we assume that a particle embedded in the fluid has a finite 
mass. The equation of motion is in general rather complicated We assume 
a simple model as 



dv x ( dip \ dv z ( dip 

- v. 



dt \dz I ' dt V dx 



(3) 



where v x = dx/dt,v z = dz/dt, and 7 _1 denotes a response time and is pro- 
portial to the mass of the particle. For the massless case, 7 becomes infin- 
ity and Eq. (2) is recovered. The volume contraction rate of the four di- 
mensional phase space (x,v x , z,v z ) is —2j, and therefore the dynamical sys- 
tem (3) is a dissipative system. We assume the same form of stream function 
ip = A/7r sin{7ra + B snx(2irb)}W(z). For simplicity, we assume W(z) — sin(7rz) 
for the convection with free boundary conditions, and W(z) = z 2 (l — z) 2 for the 
convection with rigid boundary conditions. More complicated function W(z) is 
used in the exact linear stability analysis of the Benard convection with rigid 
boundary conditions ||. 

The fixed points of the dynamical system (3) at B = are (xo, v x o, zq, v z o) = 
(1/2, 0,1/2,0), and (a*, v xl , zi, v zl ) = (0,0, 0,0) and (#2,^2,^,^2) = (0,0,1,0). 
The fixed point (xq, z$) is an elliptic point for the Hamiltonian system (2) and 
the orbit near the elliptic point is stable for small perturbation in the Hamilto- 
nian system. It is called the KAM theorem and the chaotic behaviors cannot 
be expected near the elliptic point. The fixed points (xi,Zi) and (#2)^2) are 
saddle points and the chaotic orbits appear near the saddle points for nonzero 
B in the Hamiltonian system. 

On the other hand, the fixed point (xq, v x q, zq, v z q) for Eq. (3) at B — is an 
unstable focus for finite 7 and the orbit approaches the heteroclinic orbits which 
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connects the different saddle points. Figure 1(a) displays a trajectory starting 
from (x,v x ,z,v z ) = (0.5001,0,0.5,0) for A = 1, B = 0, 7 = 20 and W(z) = 
sin7T2;. As the time evolves, all trajectories approach the saddle connections 
which connect four saddle points (0, 0), (1, 0), (1, 1) and (1,0). The numerical 
simulation was performed with the Runge-Kutta method of timestep At = 
0.0005. The orbits seem to be attracted rapidly to the saddle points. For the 




The initial condition is (x, v x , z, v z ) = (0.1,0,0.1,0). Note that the time scale 
is logarithmic. 



rigid boundary model with W(z) = z 2 (l — z) 2 , there are fixed lines at z = 
and 1. The eigenvalues around the points (x, z) = (x, 0) and (x, 1) on the fixed 
lines are calculated from 

where Sx, Sz, 5v x , Sv z are small deviations from the points on the fixed lines. 
The eigenvalues are 0,0,-7 an d — 7- It is a singular situation in a viewpoint 
of dynamical systems. The trajectories are attracted to the saddle connection, 
however, the flow near the fixed lines is very slow. Figure 1(b) displays the time 
evolution of x(t) starting from (x,v x , z,v z ) — (0.1,0,0.1,0) for A = 4n,B = 
0,7 = 20 and W(z) — z 2 (l — z) 2 . The time evolution of x(t) may seem to be 
like a limit cycle, however, the time axis is plotted with a logarithmic scale. It 
implies that it takes exponentially longer time for the particle to circulate one 
convection cell as the particle is close to the saddle connection. 

The saddle connection is unstable for general perturbations. The unsta- 
ble manifolds and stable manifolds starting from saddle points intersect and a 
strange attractor appears near the saddle points for nonzero B ||. All tra- 
jectories are attracted to the strange attractor near the saddle points and the 
chaotic diffusion is induced. Figure 2(a) is a projection of the strange attractor 
into ([x+l/2],[z+ 1/2]) plane at t = n =intcger for A = 1,B = 0.001,7 = 20 
and W(z) = sin7rz, where [y] implies the fractional part of y. The strange at- 
tractor around the saddle points (x, z) — (0, 0), (1, 0), (0, 1) and (1, 1) is shifted 
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Figure 2: (a) Strange attractor in the ([x+ 1/2], [z + 1/2]) plane for A = 1, B — 
0.001,7 = 20 and W(z) = sin7rz. (b) Stationary distribution P(z) of z for the 
same parameter values as (a).(c) Time evolution of x(t) for the same parameters. 



around (1/2, 1/2) to clarify the attractor. We have performed several numeri- 
cal simulations using different initial conditions, however, we have obtained the 
same strange attractor. The stationary distribution of the particle position is 
uniform for the Hamiltonian system (2) which corresponds to the massless case, 
however, the stationary distribution of the particle position is not uniform in the 
dissipative system (3). Figure 2(b) displays a double logarithmic plot of the dis- 
tribution P(z) of z near z = for A = 1,B = 0.001,7= 20 and W(z) = sinirz. 
The particle is attracted to the strange attractor near the saddle points and the 
stationary distribution seems to obey a power law with exponent about 0.93. 
Figure 2(c) displays the time evolution of x(t) for the same parameters from the 
initial condition (x,v x , z,v z ) = (0.5001,0,0.5,0). Chaotic diffusion is induced 
even for the very small perturbation. 

There are many "windows" of stable limit cycle attractors when the param- 
eter B is changed. If a limit cycle which runs across convection cells exists 
and it is stable, x(t) increases or decreases with a constant average velocity. 
The limit cycle corresponds to the accelerator mode in the Hamiltonian system 
(2) 0. As B is changed, the stable limit cycle disappears at a critical point. 
However, the velocity correlation is very long near the bifurcation point of the 
limit cycle. The velocity v x keeps positive or negative values for a long time 
and the direction of motion changes intermittently. Figure 3 displays such a 
constantly increasing time evolution of x(t) at B = 0.002 and an intermittent 
time evolution of x(t) at B = 0.002165 for A = 1,7 = 20 and W(z) = sin7rz. 
The diffusion constant D = lim r _ >00 ((a;(t + r) — x(t)) 2 )/2r is enhanced near the 
bifurcation point. Figure 3(b) displays the diffusion constant as a function of 
B. The average (• • •} to evaluate the diffusion constant is numerically calculated 
as double averages of a long time average and the average with respect to 100 
different initial conditions. The diffusion constant increases as 1/(B — B c ) a , 
where B c ~ 0.00216 and a = 0.5 ~ 0.6. 

The chaotic diffusion occurs both for free boundary models and rigid bound- 
ary models. The difference between the two boundary conditions is clearly 
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Figure 3: (a) Time evolution x(t) at B = 0.002 (solid line) and 0.002165 (dashed 
line) for A = 1,7 = 20 and W(.z) = sin7T2;. (b) Diffusion constants for several 
B's near the critical point of the limit cycle. 




seen, when the perturbation amplitude B is very small. Figure 4 displays the 
averge period for a particle to circulate one convection cell as a function of B. 
The period is numerically obtained as the average interval between neighbor- 
ing two times where z(t) = 0.5 and dz/dt > 0. Figure 4(a) shows the results 
for A = 1,7 = 20, W(z) — sin7rz. The average period T is approximately 
T ~ — 4.1 log 10 _B. The logarithmic dependence is characteristic of the saddle 
connection. The period increases very slowly as B is decreased. Figure 4(b) 
shows the results for A = 47r,7 = 20 and W(z) — z 2 (l — z) 2 . The average 
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Figure 4: (a) Semi-logarithmic plot of the average circulation period for W(z) = 
sin7rz. (b) Double- logarithmic plot of the average circulation period for W(z) = 
z 2 (l-z) 2 . 



period increases approximately by a power law: T ~ B -1 / 2 as B is decreased 
for the rigid boundary model. It may be due to the existence of the fixed lines 
for the rigid boundary conditions. The average period T plays a role of a unit 
time for the chaotic diffusion. For example, if the direction of motion is as- 
sumed to be randomly changed after each circulation of convection cells, the 
total steps N of random changes of direction during the time interval r is esti- 
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mated as N ~ t/T, therefore, the diffusion constant is estimated as Z? oc 1/T. 
Figure 5(a) displays numerically obtained diffusion constants at several -B's for 
A = 1,7 = 20, TT(z) = sin7rz. The diffusion constants depend on B in a com- 
plicated manner. It is partly due to the "windows" structures by the stable 
limit cycles. For example, the diffusion constant D is zero at B = 3 x 10~ 5 , 
where there exists a stable limit cycle which circulates in a single convection 
cell. However, the diffusion constant has a tendency to decrease very slowly as 
B is decreased. Figure 5(b) displays numerically obtained diffusion constants at 
several B's for A = 47r,7 = 20 and W(z) = z 2 (l - z) 2 . The diffusion constant 
decreases with a power law: D ~ B 1 ! 2 as B is decreased. The numerical results 
are consistent with the above rough argument. 




(b) 
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Figure 5: (a) Semi-logarithmic plot of the diffusion constants for W(z) — sin7rz. 
(b) Double-logarithmic plot of the diffusion constants for W(z) = z 2 (l — z) 2 . 



To summarize, we have numerically studied the chaotic diffusion of parti- 
cles with finite mass in oscillating convection. The diffusion is induced by the 
strange attractor near the saddle points. In a sence, the chaotic diffusion is 
enhanced in such dissipative systems, since there exist no stable oribits near 
the center of the convection cell and the orbits are attracted toward the saddle 
points in contrast to the Hamiltonian system. The strange attractor exists even 
for very small perturbations, since the attractor of the nonperturbed system 
is a singular one, that is, the saddle-connection. Special enhancement of the 
diffusion constant occurs near the bifurcation points of limit cycles correspond- 
ing to the accelerator modes in the Hamiltonian system. The behaviors of the 
chaotic diffusion at very weak B depend on the boundary conditions. For the 
rigid boundary conditions, the motion of the particle becomes very slow near 
the top and bottom boundaries, since the flow velocity is zero at the boundaries. 
Therefore, the average period circulating around one cell becomes very large for 
small B and the diffusion constant decreases more rapidly than the case of the 
free boundary conditions. However, it is left to future study to explain the B 
dependence theoretically. 

We would like to thank Professor S.Kai, Dr.Y.Hidaka and Mr.K.Tamura for 
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